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Lower hybrid oscillations in multicomponent space 
plasmas subjected to ion cyclotron waves 

G. V. Khazanov, 1 E. N. Krivorutsky, 2 T. E. Moore, 1 M. W. Liemohn, 3 
and J. L. Horwitz 2 

Abstract. It is found that in multicomponent plasmas subjected to Alfv6n or fast 
magnetosonic waves, such as are observed in regions of the outer plasmasphere and ring current- 
plasmapause overlap, lower hybrid oscillations are generated. The addition of a minor heavy ion 
component to a proton-electron plasma significantly lowers the low-frequency electric wave 
amplitude needed for lower hybrid wave excitation. It is found that the lower hybrid wave energy 
density level is determined by the nonlinear process of induced scattering by ions and electrons; 
hydrogen ions in the region of resonant velocities are accelerated; and nonresonant particles are 
weakly heated due to the induced scattering. For a given example, the light resonant ions have an 
energy gain factor of 20, leading to the development of a high-energy tail in the H + distribution 
function due to low-frequency waves. 


1. Introduction 

1.1. Low-Frequency Wave-Plasma Interaction 

Low-frequency Alfv6n, fast magnetosonic, and ion 
cyclotron (IC) waves are commonly found in space plasmas. 
One example is IC waves generated by ion temperature 
anisotropies in the ring current [Cornwall, 1965, 1966; 
Kennel and Petschek, 1966; Liemohn , 1967], IC waves have 
also been observed at geostationary orbit [Young et al, 1981; 
Mauk, 1982; and references therein], in the outer magneto- 
sphere [Anderson et al, 1992], at midaltitude high latitudes 
[Erlandson et al, 1990], and at ionospheric altitudes [Iyemori 
and Hayashi, 1989]. 

Low-frequency waves (LFWs) in the IC range can effectively 
interact with different components of the magnetospheric 
plasma through heating and scattering (e.g., reviews by 
Gendrin [1981, 1983, 1985] and Chang and Andre [1993]. 
Such LFWs are believed to be responsible for precipitation of 
ring current protons [Kennel and Petschek, 1966; Cornwall et 
al, 1970], formation of SAR-arcs [Cornwall et al, 1971] and 
acceleration of heavy ions [Omura et al, 1985; Tanaka, 
1985]. Noncollisional damping of IC waves may lead to the 
heating of cold electrons and ions in the outer plasmasphere 
[Galeev, 1975; Gorbachev et al, 1988; Konikov et al, 1989], 
deduced from the observed increase of plasma temperatures in 
this region [Decreau et al, 1982; Gringauz , 1983, 1985; 
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Comfort et al, 1985]. The transport of the dissipating wave 
energy from the plasmasphere along the geomagnetic field 
lines into the ionosphere below could account for the observed 
subauroral rise in the electron temperature [eg., Kozyra et al, 
1986, and references therein]. 

By interacting with charged particles, the waves influence 
the behavior of the plasma as a whole, and therefore wave 
effects must be included in corresponding models. Ganguli and 
Palmadesso [1987, 1988], Singh [1988], Singh and Torr 
[1990], Brown et al [ 1991, 1995], and Lin et al [1992, 1994] 
took into account effects from the low-frequency (LF) electro- 
static turbulence, and they demonstrated that wave-particle 
interactions lead to significant effects on the evolution of the 
core plasma distribution functions. 

In developing a mathematical model to describe plasma 
transport in the magnetosphere-ionosphere system that 
accounts for the active wave processes which occur there, we 
must develop a general scheme to include an analysis of the 
dispersion characteristics of the medium in order to choose 
suitable wave modes, a wave-particle interaction mechanism, 
and a system of hydrodynamical equations governing macro- 
scopic plasma parameters which properly accounts for the 
presence of wave-particle interactions. One basis for this 
scheme's development has been described elsewhere 
[Khazanov and Chernov, 1988; Konikov et al, 1989; 
Chernov et al, 1990; Gamayunov et al, 1991, 1992; 

Gorbachev et al, 1992]. 

When such LF waves are propagating in a plasma, a pon- 
deromotive force is produced, which leads to significant 
effects in the plasma. Boehm et al [1990] suggested that this 
ponderomotive force can be the cause of significant density 
perturbations in the auroral zone. Allan et al [1991] and Allan 
[1992] have investigated mass transport caused by the pon- 
deromotive force of hydrodynamic waves in the middle magne- 
tosphere. It was also suggested that significant plasma ener- 
gization could occur in these regions [Allan, 1993]. Li and 
Temerin [1993] and Lee and Parks [1996] have used similar 
ideas to evaluate ion energization by large-amplitude Alfv6n 
waves (£>100 mV/m; co <C2; ) in the auroral zone. 
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It was recently demonstrated by Khazanov et al. [1996] that 
large-amplitude LFWs can generate lower hybrid waves 
(LHWs) in the auroral zone and ring current region. The ion 
energization due to the LHWs may be comparable with that 
produced by the ponderomotive force of the LFWs. The pur- 
pose of this paper is to extend the analysis of Khazanov et at. 
[1996] to multicomponent plasmas and to small LFW electric 
fields. 

1.2. LHW Excitation Due to LFWs and 
Description of the Problem 

It is well known that in the electric field of a LFW, electrons 
and ions move with different velocities. The arising ion-elec- 
tron relative velocity can lead to a current instabilities and, in 
particular, LHW excitation. Depending on the frequency of the 
LFW, the magnitude of the electric field, plasma composition, 
temperatures of different plasma species, and the external 
magnetic field, strong or weak LHW turbulence may develop. 

In the case of large-amplitude LFWs, strong lower hybrid 
(LH) turbulence could be produced [Gamayunov et al 1992; 
Khazanov et al. y 1996]. Strong LH turbulence results in 
plasma heating and acceleration [Shapiro et ai y 1993], the 
formation of ion conics [Chang and Coppi , 1981; Retterer et 
al. y 1986], and possibly the saturation of Alfv6n oscillations 
in the ring current region [Gamayunov et al., 1992]. Weak LH 
turbulence can also lead to important changes in the plasma 
state; for example, Omelchenko et al. [1994] found that the LH 
"fan" instability may accelerate ions. 

The level of energy density of LFWs which leads to excita- 
tion of LHWs in a hydrogen plasma, was discussed by 
Khazanov et al. [1996]. It can be estimated by taking into 
account that the ion-electron relative velocity needed for LHW 
excitation is of the same order of magnitude as the thermal ion 
velocity. This leads to a requirement for the LFW energy den- 
sity |£ 0 | 2 > 

g Q-> gli (1) 

8; inT o)Q(o 2 pi 

where 0) Bi and co pi are the cyclotron and plasma frequencies 
for ions, respectively; nT is the thermal plasma energy den- 
sity; and ft>o * s LFW frequency. 

Khazanov et al. [1996] assumed that the plasma consists of 
only electrons and protons and that the LFW electric field is 
strong enough not only to make the relative proton velocity 
greater than its thermal velocity, but also to create strong 
LHW turbulence. From (1), however, it can be seen that the 
LFW energy density needed for LHW excitation may be lower 
in the presence of a heavy ion component. Since the magneto- 
sphere and ionosphere have multicomponent plasmas, it can 
be supposed that LHWs can be excited in such plasmas by 
LFWs with amplitudes less than that needed for a proton-elec- 
tron plasma. 

This leads us to the problem of LFW interaction with a mul- 
ticomponent plasma due to LHW excitation. In accordance 
with this problem, we will study the following questions: 

1. How does a small mixture of heavy ions influence the 
LFW electric Field needed for LHW excitation? What are the 
growth rates of the LHWs and how do they depend on the 
heavy ion population, LFW electric field, and core plasma pa- 
rameters? What are these dependencies in the kinetic and 
hydrodynamic regimes of LHW excitation? (These are dis- 
cussed in section 2.) 


2. What is the physical mechanism that determines the 
level of LHW energy density? Does the quasilinear interaction 
of LHWs with oxygen, helium, and hydrogen ions as well as 
with electrons lead to saturation of the LHW instability, or is 
it the nonlinear-induced scattering mechanism that determines 
the LHW energy density level? What is the level of LHW 
energy? (These are discussed in sections 3.1 and 3.2.) 

3. What is the influence of weak LHW turbulence on differ- 
ent plasma species? (This is discussed in section 3.3.) 

In section 4 the suppositions used and results obtained are 
summarized. Fulfillment of these suppositions and the possi- 
bility of LFW interaction with magnetospheric plasma due to 
LHW excitation are analyzed, and an estimation of particle 
energization for a magnetospheric plasma is conducted. 
Section 5 contains a short summary of the results. 


2. Linear Analysis 

2.1. Ion Drift in a LFW Electric Field 


Let us consider a multicomponent plasma containing elec- 
trons and oxygen, helium, and hydrogen ions subjected to 
LFW. In a variable electric field E (for instance, the electric 
field of an Alfv6n or fast magnetosonic wave, which has a 
component normal to the external magnetic field B 0 ) the elec- 
trons and ions will acquire different velocities, w e and w a , due 
to this field. If the frequency 0) 0 of the field is considerably 
lower than the electron cyclotron frequency, the electron and 
ion velocities will be given by [e.g., Akhiezer , 1975] 


w. 



E x Bq 


, g(E-b)b 

m e O)Q 


W “ = / H 1~T rr[fl>a a Exb + /a) 0 (Exb)xb]> (2) 

[m a [(Di a -(o 5) j 

E = E 0 e~ iav 


where b is a unit vector in the direction of the external mag- 
netic field Bq and co Ba is the ion cyclotron frequency of the 
ion species a. In a LFW the parallel electric field component 
can be small and the perpendicular drift velocities of the elec- 
trons and ions may be comparable to or even larger than the 
parallel electron velocity. This study will be restricted to the 
case where the parallel electric field and the parallel electron 
velocity can be neglected. 

In a coordinate system where the magnetic field B 0 is paral- 
lel to the z axis and the LF electric field E is parallel to the x 
axis, the relative velocity of the ion species a with respect to 
the electrons, u a = w a — w e , i s 


eE 0 w 0 


™ a (a>o -Ct>la) 


. (On 

i sin (Of)t + j — —cos o) Q t 

®Ba 


(3) 


For LFW with left-hand and right-hand circular polarization 
propagating along the z axis, the relative velocity is given by 
the expression 

u„=-7= — eE o a o [isinci) 0 / + jcos<u 0 j] (4) 

V2 m a {±(o 0 -(o Ba )(o Ba l 

If (O 0 is not near the cyclotron frequency of the light or heavy 
ions, co B < Q) 0 < 0) B , then the ratio of the velocity of the 
heavy ions u H to the velocity of light ions u L is of the order 
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u L co 0 


for waves with linear or circular polarization. For co 0 ~ co B , 
again for both cases, we get 


u 

u 


H 

T 




V 1 

) 


It may be concluded from the estimations of (3) and (4) that 
the relative ion velocity in the electric field of LFW with fixed 
frequency increases with ion mass. The reason is that the rela- 
tive velocity of the various ions is caused by the different iner- 
tia of ions and electrons in the oscillating electric field. For a 
LFW frequency near co B the resonance effect also becomes 
important. The role of the resonance denominator in (3) and 
(4) can be significant in some cases. For example, IC waves 
with a frequency near the helium gyrofrequency, generated in 
the ring current region and propagating in the direction of 
stronger magnetic field, may move to a region where co 0 is 
closer to the oxygen gyrofrequency. For waves observed by 
the Active Magnetospheric Particle Tracer Explorers (AMPTE) 
Ion Release Module (IRM) satellite [LaBelle and Treumann , 
1992], the relation (Oq/(O b was 1.2 and the relative veloc- 
ity of oxygen ions was 20 times greater than that of hydrogen. 
It follows from the above argument that in a multicomponent 
plasma, LFWs due to inertial drift of the heavy ions form a cur- 
rent across the magnetic field. 


2.2. Lower Hybrid Wave Generation 

Let us consider the LHW generation in the linear approxi- 
mation. The dispersion equation in this case has the form 
[e.g., Akhiezer , 1975] 


1 + 7T^-[ 1 + / o( JC ) e X Ze Z { z e)] + Y< 7TT [ 1 + z « Z (At)] = 0 

* v 7> a * V T a 


*l v 7> 


Z* = 


V2 k z 


Zn = 


v 7> 


<u-k x u a 


(5) 
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ne 2 2 n a el 
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m e e Q r m a e 0 



where / 0 (jc) is the modified Bessel function; Z(z a ) is the 
plasma dispersion function; and co and k are the frequency and 
the wave vector of the LH oscillations, respectively. 

In obtaining (5), it was assumed that co ^ « co 2 « co\ t 
and the motion of all ion species is unmagnetized, i.e., 
» 1 . This inequality is easily satisfied, taking 
into account that the ions with velocity u mainly generate 
LHWs with wave vector k ~ co/u, where co ~ (Q)g e co gH + )^ 2 and 
u > v Ta . The dispersion equation (5) is valid if during a time of 
order y ] (where y is the growth rate of LHW), the ion motion 
can be considered regular. This leads to the condition 


Y>0) o 


(6) 


We will restrict our analysis of (5) to the case 

c O + <<C He + <<C H + ’ where c a ~ n a / n • 

2.2.1. LHW generation in the kinetic limit. Let 

us examine the kinetic limit of LHW generation. Considering 


z e , ^h + ' z He + >:> 1’ neglecting the first term and w H+f « He+ 
in (5), and taking into account only the oxygen term in the 
imaginary part of (5), the dispersion law becomes 


co ‘ 




c H .+ 


He* 


H* 


cos 2 e+r(e) 


*V 

V 7Yt 
*>LH 


(7) 


where k\\\ e jco 2 Be <1, co LH = yjco Be co BH + , 6 is the angle be- 
tween the wave vector k and the magnetic field, and 
2.5<r(0)<4; r(0)~4. The function r(9) depends on the 
different ion species' concentrations c a , and the minimum 
r(9) is at 6, satisfying the condition m, / m e cos 2 # - c H+ . 

The growth rate of LHW is 


r=- 




CO' 


2 m * 2 v* h . 


T ' 2 

-rl 0 (x) e - x z e e~ Zr + £c a z a e' 
l e a 


( 8 ) 


where R{9) is a slowly changing function of k and 9 , 


*(*) = 


3 l m H + 9 n 

1 + -=--* — - — COS 2 9 + 
z: m. 


1 + — 

V V ) 


It is clear from (8) that a positive y can be obtained for weak 
LFW electric fields due to the oxygen ion drift (which is the 
largest drift) when in the oxygen term, co -k ± u Q + <0 and 
cu-k 1 u 0+ /k ± v Jo+ < 1 . For a small oxygen concentration, 
y> 0 in a narrow region of wave velocities, cof k^_ ~ w Q+ , with 
« 0+ > v j„+ * v j He + * * n ^is case » oxygen term is at a maxi- 
mum and the number of hydrogen and helium ions in resonance 
with LHW is small. Also, electron Cherenkov damping limits 
the generation, so LHW excitation is possible only for 

cos 2 e<\ 

T, m e 

Figure I illustrates these features of LHW generation, showing 
y >0 only for large 9 and a certain range of cof k . 

For a purely hydrogen plasma, y>0 if tu-kj_u H+ <0 and 
co- k x u H+ /£j_v T h+ < 1 , which means w H+ ~ v T h+ . This condi- 
tion leads to much greater required LFW electric fields than for 
« 0 + >v 7 y>- From (8) the maximum growth rate (neglecting 
the damping) is 



Figure 1. The normalized growth rate of lower hybrid waves 
(LHWs) as a function of the wave normal angle 
[T { m e I (T e m i cos 2 0)] and wave velocity [co/{4lk \^ )]. 
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/max 


0. 3 C Q+ 



(9) 


The restriction given in (6) leads to an inequality that is easy 
satisfied when u Q+ equals several times and c Q+ >10“ . 

Note that LHW excitation is possible if the oxygen veloc- 
ity is of the order of the thermal hydrogen velocity. Taking 
z 0+ >>1, the dispersion equation (5) can be presented in the 
form 


1 + Se ' (k, <») + <5e H * (k. co) + Se He ' (k, co) — ^ : f =0(10) 


co - k 


iv) 


where Se e , Se H+ , and 8e Hc+ are the electron, hydrogen, and 
helium dielectric permittivities, respectively. 

Assuming co = kj_u Q+ + 77 , 77 « we can find the growth 
rate from (10) 


77= . ' + " UU 

y 1 + 8e e (k, k L u 0 + ) + <5e H+ (k, k ± u 0+ ) + Se"** (k, k ± u 0+ ) 

If « Q+ « v^+ , the term <fe H+ has an imaginary part caused by 
the presence of the resonant hydrogen ions, and so one root of 
(11) leads to an instability. For co 2 pH + j k 2 v^ + »1 the esti- 
mation of the growth rate is 


Y~ 



(12) 


2.2.2. LHW generation in the hydrodynamic 
limit. Let us examine the hydrodynamic region of LHW 
generation, z a >z e » 1, x « 1 . For angles cos 2 0-m e /m 1 , 
taking into account only the relative motion of oxygen ions, 
(5) can be transformed into 



Assuming co - kj_u 0 + + 77 , where (u*=kxu 0 > is the solution 
of (13) for c 0+ = 0, the maximum growth rate of LHW comes 
out of this expr ession as 


r}= yr 


k l U Q £ 
4 


A 

o + 


(14) 


The hydrodynamic consideration is valid if » 1 , i. e.. 


w O + 


v 


7o + 



(15) 


and l77l>(W 0 (see (6)). Both of these conditions are easily sat- 
isfied. 

Using the estimations for u Q+ and « H+ (see (3), (4), and 
below), (15) may be rewritten as 


“h* >; , m H+ C ~K 1 

v Th* m O' ° i 


Note that the obtained results are valid for plane polarized 
as well as circularly polarized LFW. In the case of circular 
polarization the oxygen ions act as a beam of oscillators, with 
a rotation frequency co 0 (see (4)). The radius of rotation, 
u 0 + /ft * 0 » i s much greater than the LH wavelength, because 


u 2..2 * 2 v 2 




V 


wl 




»1 


and for this linear analysis of LHW generation the rotation can 
be neglected. 

Comparing the growth rates of (9), (12), and (14), it can be 
found that the strongest dependence on the oxygen ion density 
as well as the LFW electric field takes place in the kinetic 
regime of LHW generation in (9). Of course, in this regime the 
LFW electric field needed for instability excitation is smaller 
than in the hydrodynamic case (15). From the above deriva- 
tion it follows that in a plasma subjected to LFW the heavy 
plasma species acts as a beam of particles with a relative drift 
velocity. This heavy ion beam determines the character of the 
LHW excitation, growth rate, and restrictions on the excita- 
tion. It also follows that in a plasma with some oxygen ions, 
the LHW instability can be excited by LFW with an electric 
field amplitude less than in a purely hydrogen plasma. 


3. Quasi-linear Approximation 

The LHW instability excited by the oxygen ion drift in the 
LFW electric field leads to the development of LH turbulence. 
The growth of the LHW can be restricted owing to their quasi- 
linear interaction with different plasma species. If such inter- 
action does not lead to the instability saturation, the level of 
LH turbulence is determined by nonlinear processes, in partic- 
ular, by induced scattering of LHW by ions and electrons. 
Below, these processes are studied and the quasi-steady-state 
level of the LHW energy density is found. 

Let us summarize the characteristics of LHWs excited in a 
multicomponent plasma subjected to plane polarized LFW, 
with (Oo»0)bo+ and some other results of the previous anal- 
ysis. From (3) it follows that the main component of the rela- 
tive velocity, uo+, is directed along the y axis in the same 
phase as the LFW electric field. We will restrict our considera- 
tion to weak LFW electric fields (small E 0 ). It follows from (8) 
that in this case the LHW may be generated in a narrow region 
of angles near u 0 + (the LFW electric field needed for excitation 
of LHW is at a minimum when kjjlu 0 +)- It was found above 
that the hydrogen and helium ion motions can be treated as 
unmagnetized, so the influence of the magnetic fields on their 
motion will only lead to phase mixing. Taking into account 
that the thermal velocity of oxygen ions is much smaller than 
their relative velocity, we can consider their motion as one- 
dimensional. The region of wave velocities in which the LHW 
generation is possible is determined by the condition y>0, 
where the lower limit is given by (8) and the upper limit is 
u Q+ and is narrow for small LFW electric fields (Figure 1). 

3.1. Quasi-linear Equations 


For a purely hydrogen plasma the LHW hydrodynamic insta- 
bility is possible if « H + / v 7h + >:>1 * 11 f°U° ws f rom tJie last 
two inequalities that for a multicomponent plasma the LHW 
hydrodynamic instability occurs for a smaller drift velocity 
and LFW electric field than for a purely hydrogen plasma. 


Let us now obtain the quasilinear equations describing the 
behavior of the LHW spectral energy density and distribution 
functions of different species of plasma particles. 

The kinetic equation for oxygen ions in the presence of a 
LFW electric field is 
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^l + v .^L + i ^ cos(Uo ,%l + ^E%. 

at or m (j . av' m () . av' 


df n . 

+ “>BO*W^T = 0 


(16) 


Here / 0 + is the oxygen ion distribution function and E is the 
electric field of LHW. Introducing the relative velocity 
v = V- jE 0 cosa>Qt / B 0 , (16) becomes 


dt 


E ^- 

80 d\ m Q+ d\ 


- u 


By averaging (17) over the characteristic time of the LH 
oscillation, the quasilinear equation for the oxygen ion dis- 
tribution function can be found 


dimensional, the small component of the wave vector along 
the magnetic field k z must be prescribed. 

The changes in the LHW energy density spectrum are 
described by an equation which takes into account the energy 
generation due to oxygen ion drift and energy damping due to 
the absorption by all other species [Shapiro and Shevchenko , 
1988]: 


^■^^*1’- i-N 2 


vi^O >/* 


Oft a 

N 2 

(9vx . 

\Jk 2 v±-a) 2 


( 22 ) 


+ NLT 


The last term, NLT y takes into account nonlinear interaction 
processes of LHW with particles. In this case, it is induced 
scattering by ions and electrons [Musher et ai, 1978]. The 
system of equations (19)-(22) describes the evolution of dis- 
tribution functions and LHW electric field energy density. 


where u = eEo/ma>o , v=Vy, and b(x) is a delta function. 
Equation (18) follows from (17) if the oxygen ions drift along 
the x axis and the magnetic field terms are neglected. The first 
assumption is valid with an accuracy of 0) bo +/(Oq« 1, the 
second with an accuracy of v To +/w«l. Supposing the LHW 
spectrum to be one-dimensional, evaluation of the integral in 
the right-hand side of (18) yields 





# 0 + _ 7te 2 
* " 


j l £ *| 2 fro* 

d\ V <9v 


(19) 


where |E*| 2 is the LHW spectral energy density. It was 
assumed in (19) that |v- dcojdk \ , arising after the integration, 
is equal to v. 

Taking into account the phase mixing and supposing the 
LHW spectrum is one-dimensional, the quasilinear equations 
for the resonant part of the hydrogen and helium distribution 
functions [Shapiro and Shevchenko , 1988] can be described in 
the form 


3,2 LHW Energy Density 

Our purpose here is to determine the LHW energy density 
and the influence of LHW on oxygen, hydrogen, and helium 
ions and electrons. We will avoid the beginning stage of LHW 
instability development and restrict the analysis to a time 
long enough for significant changes to occur in the initial dis- 
tribution functions of the particles. 

Let us examine (20) describing the distribution functions of 
hydrogen and helium ions. Owing to the interaction with the 
resonant LHW, these ions are accelerated, and, for large times, 
it can be supposed that the resonant ions have velocities such 
that vx»ft)/fc_|_. Then (20), for (o/k^- 0, have an asymptotic 
solution for t— that does not depend on the initial distribu- 
tion function and is self-similar in form [Shapiro and 
Shevchenko, 1988], 


/, 


5 ”H + (He + ) 

H 1 He+ ) = r(2/5)f 2/5 


exp 


[-vi/r 



(23) 
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j_JL f 

Vl 1 L * 2 V* vi-® 2 Vl 


^h*(hc*) 

<9vj_ 


dk (20) 


Equations (20) are integrated over the azimuth angle in the 
plane perpendicular to the magnetic field assuming axial 
symmetry of the distribution functions. It follows from (20) 
that the resonant interaction for unmagnetized particles 
(Ct)=kxVx) is possible for hydrogen and helium ions with 
vx^fiV^x- 

The resonant condition for electrons may be satisfied only 
for Cherenkov interaction. The quasilinear equation determin- 
ing the electron distribution function evolution is [Shapiro 
and Shevchenko, 1988] 


%e_ = 


Tie 


dt 


dk 


dv f 


N 2 

4{x) 

k 2 

da 

Vz ~lk 


<9v. 


(21) 


where J 0 (x) is the Bessel function. In (21), which is one- 


where T(2/5) is the Gamma function and nn + (He + ) * s the number 
of resonant particles. 

From (23) it follows that the distribution functions of 
hydrogen and helium ions, and therefore the energy absorp- 
tion, only weakly depend on fast energy density oscillations 
and thus are determined by the slow changes in the average 
energy density. For the thermal velocity of helium 

ions is half that of the hydrogen ions, so the number of helium 
ions in the region of resonant velocities (for Maxwellian 
initial distribution functions) is (rc He + / n H + )exp(-3v 2 /v^ + ) 
times less than the number of hydrogen ions. This means tnat 
the absorption by hydrogen ions is much greater than by 
helium ions. 

From (21) it follows that the electron distribution function 
is one-dimensional. Therefore, for large times, a plateau in the 
region of resonant velocities along the z axis will eventually 
form and the energy absorption will cease. Thus, for large 
times, only the term describing the energy absorption by 
hydrogen ions needs to be kept in (22). 
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Note that the maximum level of LHW energy density is 
limited by the nonlinear term in (22). This energy level can be 
estimated from the balance of the quasilinear oxygen term and 
the nonlinear term. 

The damping rates due to induced scattering by electrons and 
ions in the case under consideration are of the same order of 
magnitude [Musher et al, 1978], 

yNL _ £ ( 24 ) 

W LH ®LH nT 


where 


RjtJ 1 *' 8?r 



From the equality of the linear growth rate y (the maximum 
value of the quasilinear term) and the nonlinear damping rate 
y NL , the maximum value of the LHW energy density can be 
estimated, 


N 2 ?W.h« nT 

8?r «L©pH* 


(25) 


If we suppose that the system of LHWs and particles pro- 
gresses to a steady-state condition due to quasilinear inter- 
actions, then we can consider the slow evolution of the 
oxygen distribution function and the LHW spectral energy 
density. 

Let us investigate the solution of (20) and (22) with the 
form 


fo = F(v,t) + a(v)cosco 0 t + b(v)sin(O 0 t 
N 2 = e k (v,/)+c(v)cos£u 0 r + d(v)s\na) 0 t 


(26) 


where F and e k are slowly varying functions of time. We will 
assume that the growth rate determined by the slowly 
changing part of the oxygen ion distribution function F and 
the damping rate caused by hydrogen ions are less than (O 0 . 
This damping rate can be estimated from the hydrogen term in 
(22) using the distribution function (23) and the LHW energy 
density (25) and is in accordance with this assumption. 

In this case, the equation for F may be written in the form 

dF ne 2 u 2 d 2 e k d 2 F ne 2 d e k dF __ 

dt m 2 + 2 dv 2 v <9v 2 m dv \ d\ 


The second term on the left-hand side shows that the asymp- 
totic influence of the LFW electric Field on the oxygen ion 
motion is determine by the "friction" of the ions with the 
LHW electric field. Equation (22) for the slowly varying part 
of LHW energy density e k has the form 



7tm H+ 




o + 


O) 3 dF na> 4 r d\ ± 

k 2 n k dv k 3 n ^ v^ d\± * 


(28) 


As was discussed above, only the absorption by the hydrogen 
ions is taken into account in (28). 

We will suppose below, as is usual in a quasilinear approx- 
imation, that the LHW spectrum has a sharp maximum at the 
wave vector k~co/u. Equation (27) can then be reduced to the 
following form: 

dF | n i ^ 1 d 2 F d 1 dF 

dr 2 d \ 2 v dv 2 dv v dv 


where 


t = ^ r je k dt 

m O + o 

Let us estimate the role of the different terms on the time 
behavior of the oxygen ion distribution function. If we 
neglect the quasilinear term, then the asymptotic self-similar 
solution can be found. It does not depend on the initial condi- 
tions and can be calculated by the perturbation method as a 
series of powers of tj-v 5 rK The term which determines the 
LHW evolution for t—>oo is F~vr 2/5 . Conversely, if we neglect 
the term with « 2 + , the analogous quasilinear solution is 
F~v-V 4/3 . The first expression shows that the LFW electric 
Field causes a positive value for the term dF/dv and leads to 
LHW generation. From, the second expression for F it follows 
that the quasilinear interaction gives a negative dF/dv term. 
Such a result means that by neglecting the term with the LFW 
electric Field, we reduce the problem to a completely different 
one with the oxygen particles being accelerated by the LHW 
electric field. 

We will suppose for further estimations that, for large 
times, the behavior of the derivative, dF/dv , is determined by 
the second term in expression (29). This means that the 
energy flow to the plasma is directly related to the LFW elec- 
tric Field and leads to 


dF - V - H 

dv u w 


(30) 


Taking into account (23) and (30), the equation for the LHW 
spectral energy density (28) may be reduced to the form 


d 2 D 
dx 2 : 
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where and 


D = 


ne 


m 2 .(O 0 u 2 


X 


je k dt 

o 


and we have assumed (of k ~ u. 

In (3 1 ), the asymptotic dependence of the LHW energy den- 
sity is different for the oxygen and hydrogen terms (the First 
and second terms on the right-hand side, respectively). The 
oxygen term yields E ~ (f-/o) 3/2 , a growth of energy, while 
the hydrogen term has the form E~{t § -/) 2/3 , a damping of 
energy. Examining the powers shows that the oxygen term 
will dominate for large times and the LHW energy will grow. 

The numerical analysis of (31) gives the same results. In the 
calculations below a maximum time a) Q t - 10 3 is used, suppos- 
ing that ion velocities are 1 0 5 - 1 0 6 cm/s, the region of inter- 
action with LFW is 10 8 -10 9 cm, and co 0 is a few s' 1 . Figure 2 
shows that the solution of (31) weakly depends on the initial 
conditions. The quantity h/n can be changed in magnitude 10- 
20 times without signiFicant energy changes, allowing a wide 
range of concentrations of resonant hydrogen ions (Figure 3). 
The dependence on oxygen ion concentrations is more signif- 
icant, however. If ti 0 Jn is changed from 10' 2 to 3xl0" 2 , the 
LHW energy growth for 0) 0 t ~ 10 3 changes 20 times. Owing to 
the difference of the numerical coefFicients in the terms on the 
right-hand side of (34), the exact dependence of the LHW 
energy on r in the derivative in (31) is not very signiFicant 
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Figure 2. The normalized LHW energy density e/e 0 from (31) 
for initial conditions of e/e 0 =Q.5 (solid curve), e/e o =0.05 
(dashed curve), and e/e o =0.005 (dotted curve), as a function of 
normalized time (0 o t y with h/n- 0.005 and Z)' 2/5 . 


(for co 0 t < 10 3 ). Even for a time dependence of dFjd\ ~ t - 3/2 , 
the energy grows for (D$t - 10 3 (Figure 4) and the quasi-steady- 
state condition cannot be achieved. It follows from the above 
analysis that the quasilinear interaction does not lead to the 
saturation of the LHW instability, so the level of LHW energy 
is determined by the nonlinear process and is given by (25). 

3.3. Plasma Energization 

Now we can estimate the plasma energization produced by 
LFW due to LHW excitation. As it was found before, the quasi- 
linear interaction of LHW with helium ions is small owing to 
the large difference between the wave phase velocity and the 
thermal velocity of the particles. Since the electrons are mag- 
netized and the region of LHW phase velocities is small, the 
quasilinear interaction for electrons also can be neglected. The 
main quasilinear effect is LHW-hydrogen ion interaction. The 
LH turbulence also gives energy to the plasma due to induced 
scattering by ions and electrons. Using the spectral energy 
density from (25) and the hydrogen ion distribution function 
from (23), we can estimate the resonant hydrogen ion energy 


and the dependence of this energy on the time of interaction 
with LHW, 
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(32) 


where h is the number of resonant hydrogen ions. 

Let us estimate the energy flow Q in the plasma particles due 
to induced scattering by ions and electrons. For (0-3-46)^ 
the result is of the same order of magnitude for both species 
[Musher et ai> 1978]. For the stationary state, 

jyW k dk = Q (33) 


where Wfc is the spectral energy density of LH oscillations. 
Using (25), it is found that 


Q y 2 (o 
nT m\ e 


(34) 


4. Discussion and Results for Magnetospheric 
Plasma 

Let us discuss the main results and the assumptions that are 
used in the analyses of LFW interaction with a multicompo- 
nent plasma due to LHW excitation and the possibility of their 
applications to the magnetospheric plasma. 

In obtaining the LHW frequency (8), it was supposed that 
the LFW frequency is small compared with the LHW frequency 
and growth rate. If the LFW is an ion cyclotron wave (ICW) 
with a frequency near the helium cyclotron frequency [Young et 
al , 1981] or oxygen cyclotron frequency [Fraser et al, 1992], 
the frequency of LHW is (0 = 2g> lh (7), Y is determined by (9), 
and the condition y > f° r 58 a> fiHe + * eac * s t0 the i nec l ua l' 
ity <V w o +/ v^ »10 -2 . For the oxygen drift velocity 
w Q+ ~3 v T h+ this condition is satisfied even for an oxygen ion 
concentration c Q+ -1%. For (O 0 “ the inequality is 4 
times weaker. 



Figure 3. The normalized LHW energy density from (31) for 
the initial condition e/e o =0.005 as a function of (0 o t with 
D 2t 5 . The solid curve is for n/n=0.1, the dashed curve is 
h/n- 0.05, and the dotted curve is n/n=0.005. 



Normalized Time, cj.t 

Figure 4. The normalized LHW energy density e/e 0 from (31) 
for the initial condition e/e o =0.005 as a function of 0) o t with 
h/n- 0.005. The solid curve is for Z)‘ 7/5 , the dashed curve is 
£)" 4/5 , and the dotted curve is D’ 2/5 . 
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Another assumption used in obtaining the LHW growth rate 
is about the uniformity of the plasma, external magnetic field, 
and the electric field of LFW. The wavelength of the LHW can 
be estimated as A -2 7rw 0+ /<y LH , and even for L-7 it is of 
order A-10 5 cm. Therefore this assumption is usually fulfilled, 
because the 1CW wavelength is greater than 10 7 cm. 

Let us discuss the LFW electric fields needed for LHW excita- 
tion. LHW excitation in a plasma with an oxygen mixture of 
1-3% is possible if the oxygen ion drift velocity is 3.5-4 

times greater than the thermal velocity of hydrogen ions. In 
this case, the LHW Cherenkov damping on hydrogen ions is 
small and the development of the instability is possible. 
Taking into account that the instability develops on a 
timescale much smaller than the LFW period, the LHW genera- 
tion is possible even if the oxygen ions have such a drift 
velocity only during a short part of the LFW period. 

As was found in section 2.1 in (3), the drift velocity of a 
particle in a LFW electric field depends on the mass of the par- 
ticle and the ratio of the LFW electric field frequency to the 
particle's cyclotron frequency. Expressions (3) and (4) are cor- 
rect if \(Oo-(0 Ba \ »k 0 v Tay where k 0 is the wave vector of 
LFW. Owing to the large LFW wavelength, the right-hand side 
of this inequality is small. For ICW with a frequency of the 
order of the helium ion cyclotron frequency and an oxygen ion 
temperature 1-2 eV, &o v 7 b + ~ SxlCT 4 . Therefore (3) and (4) are 
correct for ICW frequencies close to the particle's gyrofre- 
quency. 

Now let us estimate the possibility of LHW generation 
using observed LFW magnetic and electric fields. The plasma 
composition and the hydrogen ion temperature are not given 
in the papers cited below, so it was assumed that oxygen con- 
centration is of order 1% and the hydrogen ion temperature is 
l eV. First, let us present the cases to be discussed. 

Young et al. [1981], for an event observed on board GEOS 1 
(August 16, 1977; L-7; 1200 LT), found that the electric fields 
were 2.9 mV/m (1120 UT) and 4.3 mV/m (1215 UT); the wave 
frequency was f= 0.6 Hz, the frequency range was Af~0. 1 Hz, 
and the helium gyrofrequency / He+ = 0. 45 Hz . 

The event observed by ISEE 1 and 2 (August 22, 1978; L-7, 
1700 LT) is described by Fraser et al. [1992]. For this event 
the electric field was found to be 1.4 mV/m; the wave frequency 
f= 0.12 Hz, and the oxygen gyrofrequency / Q+ =0.06 Hz. 

In the case described by LaBelle ei al. [1988] and LaBelle 
and Treumann [1992] (AMPTE/IRM observation, June 6, 
1985; L-4.5; 1547 UT) the electric field in full is not given. 
Our estimation based on this data leads to an electric field 
amplitude of -1.8 mV/m. The wave frequency in this case is 
close to the oxygen ion gyrofrequency, /=0.15 Hz, 
/ 0 - =0.125 Hz; and the frequency region is Af~Q.\5f. 

In all of these cases, the oxygen ion drift velocity is 3.5 or 
more times greater than the hydrogen ion thermal velocity, 
and so LHW can be excited. In the last two examples the reso- 
nance denominator in the expression for the oxygen ion drift 
velocity (3) is significant. Taking into account that the listed 
observations were obtained in the region near the equator for 
L- 5-7, it can be found that the difference cuo - &BO+ does not 
change significantly in a region of order R E along the mag- 
netic field lines. 

The next supposition used is that the LFW is monochro- 
matic. In two of the cases cited above ( Fraser et al. [ 1 992] and 
LaBelle et al. [1988] and LaBelle and Treumann [1992]) it is 


quasi-monochromatic with Afff~QA. Young et al. [1981] 
reported that the ICW emissions observed on board GEOS 1 
and 2 "are often quasi-monochromatic (Afff -0.1 -0.2)." 

From the previous analysis it can be concluded that LHW 
excitation due to LFW activity is a possible mechanism of 
LFW interaction with multicomponent magnetospheric 
plasma, at least at L> 5 near the equatorial region. As discussed 
by Khazanov et al. [1996], the same mechanism of interaction 
is possible in the auroral zone. 

Now we will estimate the plasma energization due to the 
discussed mechanism for magnetospheric plasma. For this cal- 
culation, let us take the magnetospheric plasma in the region 
of the plasmapause (L-4-5) with the densities of the hydrogen, 
helium, and oxygen ions in the ratio 100:10:1, respectively. 
The plasma is assumed to be isothermal with a temperature of 1 
eV. For the LFW with the frequency 0.15 Hz and oxygen ion 
gyrofrequency / o +=0.125 Hz, the needed electric field ampli- 
tude of LFW, Eq (see (8)), is 1.8 mV/m. 

The LHW frequency in that region, given by (7), is 
(0 - (0.5-1) x 10 -3 s‘\ and the growth rate is proportional to 
the oxygen ion concentration cq+ (9), y~5c 0+ <y. The wave- 
length of LH oscillation is A = 2/m/ co - 10 4 cm, and the 
envelope of generation is zlA-A. 

We will suppose that the LFW electric field as well as the 
external magnetic field and plasma parameters are the same in 
a region of order ^~10 8 -10 9 cm. Then the LHW turbulence in 
this region can be treated as uniform. The tracing time for 
hydrogen ions for such a region is of order 1 0 2 - 1 0 3 s. It can be 
taken that during this time the LFW and LHW activities are 
quasi-steady. Then, from (32), the energy obtained by the res- 
onate hydrogen ions during the tracing time can be found. It is 
of order 20-40 eV. The number of resonant hydrogen ions can 
be estimated supposing that their initial distribution is 
Maxwellian. It is equal to the number of particles with veloci- 
ties above the minimum LHW phase velocity. In the case under 
consideration, the LHW phase velocity is 3.5 v T h + and the 
number of resonant hydrogen ions is -0.5%. The energy of the 
core plasma, found from (34), is doubled for the period of 
interaction owing to induced scattering of LHW. 

5. Summary 

The preceding analysis demonstrates that LFWs in a multi- 
component plasma can generate LH oscillations. The energy 
density of the LFW electric field needed is significantly less 
than in the case of a purely hydrogen plasma. The oxygen ions 
play the role of a beam, and in the linear approximation the 
concentration dependence of the growth rate and the kinetic 
and hydrodynamic limits of LHW generation are the same as 
those of a beam instability. The quasilinear effects do not lead 
to saturation of the instability, and the energy density level is 
determined by the nonlinear processes of induced scattering. 

In the region of resonant velocities the hydrogen ions are 
accelerated and the concentration of accelerated ions weakly 
depends on the oxygen ion concentration. The nonresonant 
particles are weakly heated owing to induced scattering of LH 
oscillations. 

For a given region of the magnetosphere (L-4-5) with LFW 
electric fields of 2 mV/m, plasma temperatures of 1-2 eV, and 
an oxygen ion concentration of -1%, the resonant hydrogen 
ions (-0.5%) obtain energies of 20-40 eV. Therefore the dis- 
cussed mechanism may lead to the growth of a tail in the 
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hydrogen distribution function in the presence of low- 
frequency activity. 
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